Magnetization plateaus in weakly coupled dimer spin system 
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I study a spin system consisting of strongly coupled dimers which are in turn weakly coupled in 
a plane by zigzag interactions. The model can be viewed as the strong-coupling limit of a two- 
dimensional zigzag chain structure typical, e.g., for the (oc)-planes of KCUCI3. It is shown that the 
magnetization curve in this model has plateaus at | and | of the saturation magnetization, and an 
additional plateau at | can appear in a certain range of the model parameters; the critical fields 
are calculated perturbatively. It is argued that for the three-dimensional lattice structure of the 
KCUCI3 family the plateaus at 4 and | of the saturation can be favored in a similar way, which 
might be relevant to the recent experiments on NH4CUCI3 by Shiramura et al., J. Phys. Soc. Jpn. 
67, 1548 (1998). 
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I. INTRODUCTION 

In the last couple of years there has been a growing 
interest to the phenomenon of magnetization plateaus in 
one-dimensional (ID) spin systems .lTcj In a variety of 
models, the magnetization per site m z as a function of 
the applied magnetic field H exhibits plateaus at certain 
values of m z ; at those plateaus the magnetization in units 
of the saturation M = m z /S (here S is the spin of a 
magnetic ion) is "locked" at some rational number (at 
least, up to now no indications of the possibility to have 
plateaus at irrational M were found). Oshikawa et aln 
have shown that in purely one-dimensional systems the 
allowed values of M, at which the plateaus are possible, 
are given by the following condition 



QS(1-M) eZ, 



(1) 



where Q is the number of magnetic ions in the mag- 
netic elementary cell of the ground state (which can 
be different Appm the elementary cell prescribed by the 
Hamiltoniaroj). This condition is necessary but not suf- 
ficient; in other words, if the plateau exists, it has to 
be at one of the values of M determined by (p]), but in 
principle the plateaus do not have to exist at all values 
of M which are allowed by (|]). Validity of the result 
(|j) is supported by a considerable amount of numeri- 
cal and analytical data, as well as by the recent experi- 
mental observations of the M — | plateau in the S = 
1 bond- alternating chains realized in [Ni2(Medpt)2(/i- 
ox)(//-N 3 )]C10 4 • 0.5H 2 O. 

The present paper is motivated, by the recent high- 
field magnetization measurementsEj in the double-chain 
S = i compound NH4CUCI3 which have revealed well- 
pronounced plateaus at M = \ and M = | in the tem- 
perature range T < 1.5 K. At room temperature, the 
crystal structure of this material coincides with that of 
KCuCbj, and the main feature is presence of the dou- 
ble chains of edge-sharing CuCbo octahedra, separated 



by NH+ ions (see figs. 2,3 of Ref. |T| or fig. 1 of Ref. 
|l7| ) . Double chains composed of magnetic Cu 2+ ions can 
be viewed as spin-i bond- alternating antiferromagnetic 
(AF) zigzag chains with two different nearest-neighbor 
exchange constants J 2 , J3 and next-nearest-neighbor ex- 
change Ji (see fig. 1 of Ref. [l7|, and Fig. [I] of the present 
paper); in what follows, the word "chain" is used in the 
meaning of a zigzag (double) chain. If one assumes that 
the essential physics in NH4CUCI3 is determined by sin- 
gle chains, then the appearance of plateaus at \ and § 
of saturation means that due to a certain rather compli- 
cated spontaneous symmetry breaking the magnetic ele- 
mentary cell contains Q = 8 Cu 2+ ions instead of Q = 2 
as suggested by the Hamiltonian symmetry; still, then 
absence of plateau at M = \ which is characteristic for 
a single chairo is puzzling. 

Values of the exchange constants in NH4CUCI3, as 
well as_in_.the other materials of the so-called KCUCI3 
family,EjO are jaresently not known. There are cer- 
tain argumentsEJ based on fitting the susceptibility 
data for the isostructural compound KCUCI3 which sug- 
gest that one,-ef the nearest-neighbor couplings (J3) 
is dominatingJia On the other hand, neutron scatter- 
ing data for KCUCI3 strongly suggest that interchain 
interactions are important .E3 From the susceptibility 
datajlllll3 one can guess that the intrachain couplings in 
NH4CUCI3 are much weaker than in KCUCI3 and thus 
interchain couplings should play even more important 
role; this assumption also fits well to the fact that the 
zero-field ground state of NH4CUCI3 is magnetic at low 
temperaturelJj 

In the present paper I show that interchain couplings 
can considerably affect the number and positions of the 
plateaus. To demonstrate that, I consider a simple model 
describing a system of zigzag chains coupled in a plane; 
within the effective spinless fermion model valid in the 
limit of weakly coupled dimers it is shown that this model 
can exhibit magnetization plateaus at |, ^, and | of 
the saturation. The critical fields are calculated using 
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the lowest-order perturbation theory around the "atomic 
limit" of the model. Physics of the plateau at M = | 
is determined by single chains and is almost exactly the 
same as studied in Ref. ||, while the plateaus at M = |, | 
appear solely due to the two-dimensional interchain cou- 
plings. Generally speaking, the number and positions 
of the "additional" plateaus caused by interchain inter- 
actions strongly depend on the topology of the specific 
lattice. Possible relevance of this mechanism to the ex- 
periments on NH4CUCI3 is discussed: it is argued that 
the three-dimensional lattice structure of NH4CUCI3 nat- 
urally favors the appearance of plateaus at i and § of the 
saturation. We also discuss possible reasons of absence 
of the other plateaus in this material. 

The paper is organized as follows: in Sect. [H] the effec- 
tive model for weakly coupled dimer system is introduced 
and the magnetization curve of the two-dimensional 
system of coupled chains is studied, Sect. Ill is de- 



ft = flV'ly + U V I',, II, | I., + t) (cl i C i+ l t j + h.C.) 



voted to the analysis of characteristic features of three- 
dimensional coupled chain system typical for NH4CUCI3, 
and Sect. IV contains discussion and concluding remarks. 



II. TWO-DIMENSIONAL COUPLED CHAIN 
MODEL 



ij ij 



Z-^i v y 

ij 



ij 

t' X' ! ''</''•••'• 1 + 4j c i-l,j+l + h - c -) ■ 



(3) 



Here = cj-Cy, and n,j = 1 corresponds to the \t + ) 
state on the (ij)-th dimer, while absence of particle 
(riij — 0) corresponds to the singlet. The magnetization 
per site in the saturation units M coincides now with 
the concentration of fermions. The effective interaction 
constants U , U' , hopping amplitudes t, t' , and chemical 
potential fj, are given by 

U = (2Jj + J 2 )/4, t = (2Jj - J 2 )/4 

U' = J'/A, t' = -J'/4, v = J 3 -gti B H. (4) 

It is useful to note the particle-hole symmetry of the ef- 
fective Hamiltonian: written in terms of hole operators, 
(^) preserves its form (up to a constant shift), with the 
change only in the value of the chemical potential 



M <-> Mhoio = -fi-2U - 4U' . 



(5) 



Consider a system of zigzag chains coupled in a plane 

r frustrating zigzag interaction J', as shown in Fig. 
Such a structure is typical, e.g., for (ac)-planes of 
KCUCI3. Let us further assume that the exchange con- 
stant J3 is antiferromagnetic (AF) and much stronger 
than all the others, then the system can be viewed as a 
two-dimensional arrangement of dimers weakly coupled 
with each other. Assume further that the external mag- 
netic field H is applied in the z-direction. The model is 
described by the Hamiltonian 



H 



j) a(i,j) 



Jl(s («) . s f+hj) 

j,( S m . s (id+i) . 



J 2 s\ 



g(i,j) , J + l) < 



g(j, B H ■ (S\ 



(i,3) 



: (i,3) 



(2) 



where g is the Lande factor and /is is the Bohr magneton. 
If one neglects weak interdimer interactions J±, J2, J' 
completely, than the lowest states at each dimer are the 
singlet \s) and the S z — +1 triplet \t + ). We will be 
interested in the regime of strong fields g^sH ^ J3, then 
the other two triplet states |to), \t-) are high in energy 
and their contribution can be neglected, in spirit of Refs. 
|9|,|l2],[l4|. In this way one gets a reduced Hilbert space with 
only two degrees of freedom per dimer, and the problem 
can be formulated in terms of spinless fermions living on 
the asymmetric triangular lattice as shown in Fig. |^, with 
the following effective Hamiltonian Ti: 



In what follows we assume that the effective interaction 
constants U, U' are positive, i.e. the interaction is repul- 
sive. 



A. The "atomic limit" 

Consider first the "atomic limit" of the model (||), i.e. 
put hopping t, t' to zero. Then the problem becomes 
equivalent to the generalized Ising model in magnetic 
field, a class of models which is well studied in the theory 
of alloys (see, e.g., Ref. ^l] and references therein). Gen- 
erally, in this limit the interaction energy per site Ei nt 
is a piecewise linear function of the concentration M, 
and the curve E- lnt (M) has "kinks" (jumps of the first 
derivative) at certain values of M. When the magnetic 
field increases, the chemical potential fi becomes nega- 
tive, and particles start being "pumped" into the system 
until this process is stopped by the interaction. Kinks 
in _E; nt (M) correspond to plateaus in the magnetization 
curve M{H), and the width of a plateau is proportional 
to the magnitude of jump in the first derivative. The 
atomic limit determines the possible number and posi- 
tions of plateaus; hopping generally tends to smear out 
the plateaus, and below we will study this effect pertur- 
batively. 

In the atomic limit the model (||) is equivalent to the 
so-called anisotropic triangular nearest-neighbor Ising 
model in. external field, and its ground state is known 
exactly.t3 Up to M = | one can insert particles into the 
system without loosing any interaction energy (see Fig. 
Thus, 
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E- mt (M) = 0, 



M < - . 
~ 3 



(6a) 



Using the particle-hole symmetry, one immediately ob- 
tains from the above the behavior of E; n t above M = 



3 - 



E int (M) = (2U + AU')(M- -) 



M >_l 



(6b) 



< 



In the "intermediate" interval of concentrations 
M < | the situation depends on the ratio of interac- 
tion constants U and U' . If U' < U, then the lowest 
energy cost for inserting new particles above M = -| is 
3U' per particle, which is determined by the multiparticle 
process of inserting at once a row of particles as shown in 
Fig. ^a. This process leads from the state with M = | to 
the state with M = i , and another multiparticle process 

], taking 



leads from M = \ to M 



shown in Fig. 

the energy of 2U + U' per particle. 

For U' > U, similar multiparticle processes exist which 
govern the increase of concentration above M = i and 
M = i, with the energy cost of 3U and AU' — U per 
particle, respectively (one can obtain the corresponding 
pictures by rotating the particle configurations of Fig. 
||a,b by 7r/3 but keeping the solid and dotted lines fixed). 

Taking all that into account, one obtains the following 
result for the interaction energy as a function of M in 
the [^, |] interval: 



f 3E/'(M-i), 

\ (2U + U')M - U, 



i < M < \ 
\<M<\ 



(6c) 



for U' <U, and 

_f3tf(M-§) 
hmt ~ \ (AU 



U)(M 



- < M < - 
-) — — 1<M<1 

2 1 2 ' 2 3 



for U' > U. The resulting behavior of E- mt {M) is 
schematically shown in Fig. [|. Thus, generally three 
plateaus are possible, at M ==h, \, and |. In the limit 
of zero interchain interaction U = only the plateau at 
M = i survives, so that one may think of this plateau 
as coming from single chains, while plateaus at ^ and 
| appear due to 2D interchain coupling. The plateau at 
M = \ disappears when 2D interactions become isotropic 
(i.e., U' = U). 

It should be remarked that the positions of those "ad- 
ditional" plateaus strongly depend on the topology of 
the lattice. For instance, if one considers a trellis lat- 
tice shown in Fig. |^ (this type of interchain coupling is 
realized, e.g., in Sr n _iCu n +i02 n familja), then spinless 
fermions live not on a triangular, but on a square effec- 
tive lattice, and the only possible plateau would be at 
M 



2 ' 



B. The effect of hopping 

Since in absence of hopping we know the solution ex- 
actly, it is natural to try to take hopping into account 



perturbatively. Then one can determine the critical fields 
marking the beginning and the end of each plateau. 

To justify the use of perturbative approach with hop- 
ping term as a perturbation, one has to demand that 
J2 is close to 2Ji, i.e. that a single zigzag chain is close 
to the Shastry-Sutherland line,E£l and that the interchain 
coupling is much smaller than J\\ 



|2Ji - J 2 | < Ji, J'<Ji 



(7) 



Comparing ([?]) with (|4|), one can see that it does not 
make sense to consider the case U' > U, and in what 
follows we will assume that U' <C U. 

Let us start with the M = ^ structure shown in Fig. ||, 
where one third of the total number of sites N is occupied 
by particles. Increasing of the magnetization from M = 
^ is determined by the multiparticle process of inserting 
a row of L new particles, L 3> 1 (see Fig. |a). Then 
there are no first-order corrections to -E/v/3+L, and in the 
second order the gap A C 3 = {E N / 3+L — E N / 3 )/ L reads as 



A. 



u'<u 

c3 



H + 3U' + 

2t 2 - I2t' [ 
U + U' 



2U-W 



At 2 

IT 



u> 



(8) 



Zero of the gap A C 3 gives the critical field H C 3 . 

Similarly, for U' < U, increase of the magnetization 
from the M = ^ state is determined by the multipar- 
ticle process shown in Fig. |3f3, which, up to the sec- 
ond order, yields the following expression for the gap 
A C 5 = {E N /2+L - E N / 2 )/L: 

A c5 = /i + 2U + U' + {At 12 - M 2 )/U + At' 2 /{U + U') 
+ 3t 2 /U' -Qt' 2 /{2U -U') + 0{l/L). (9) 

To obtain the critical field H C 2, one has to consider 
the set of degenerate states with one particle taken out 
of the structure shown in Fig. 0; the degeneracy is lifted 
only in the second order, giving the following dispersion 
£c2(&) = E N / 3 „i(k) — E N / 3 of the N/3 — 1 excitation 

g 

£ C 2(fc) = -H+ jjj[i 2 + t' 2 + tt'(cosk 1 +cosfc 2 )] (10) 
8t' 2 r „ , N1 3t 2 I2t' 2 

and the gap A C 2 = mine C 2(fe) is given by 

A c2 = -ix - M 2 /U' - I2t' 2 /{U + U') 

+ 8(|i| - \t'\) 2 /U' + m l2 /U, \t/t'\ > 2U'/U, 

A c2 = -fi + M 2 /U' - I2t' 2 /{U + U') (11) 
- At 2 U/U' 2 + 8t' 2 /U', \t/t'\ < 2U'/U . 

Finally, it is easy to calculate the value of the first criti- 
cal field H r \ which-indicates the point where the injection 
of triplets starts.EJ In the vicinity of H c \ the density of 
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triplets is very low, so that one can just neglect the in- 
teraction terms in (|^) . The hopping part of the effective 
Hamiltonian yields the dispersion E\ — Eq = e(k), 

e(k) = /i + 2t cos k\ + 2i'{cos + cos(fci — fca)} • 

The gap A c i = mine(fc) is given by 



H-2t- (t') 2 /t at t > 
fj, + 2t-4\t'\ at t < 



(12) 



All the other critical fields now can be obtained by 
exploiting the particle- hole symmetry (|^), which connects 
the following pairs: 



Mc3 <->■ ^c6, 



A*c5 ' 
Mc2 



Mc4, 
' Mc7 • 



Zeros of the gaps given by the formulas (g), (||), ([TT]) , 
( |l2] ) , together with the definition of /i as stated in (Q) , de- 
termine the critical fields. It should be remarked that in 
the derivation of ( [l2] ) we did not use perturbation theory, 
so that the expression for H c \ is exact within the effec- 
tive model (||). It can be improved, however, by including 
the processes involving the other two triplet states which 
were discarded in our consideration, in spirit of the de- 
generate perturbation theory (see, e.g., Ref. |25|). 



III. THREE-DIMENSIONAL NH 4 CuCl 3 -TYPE 
LATTICE 

Now let us proceed to a more complicated model, 
which could be useful for understanding the physics of 
magnetization process of NH4CUCI3. Consider a three- 
dimensional model which can be imagined as set of two- 
dimensional systems ( "layers" ) from Section [n] stacked 
on top of each other, as shown schematically in Fig. |t]. 
Besides exchange interactions Ji.2,3, J' inside a single 
layer, we have now two additional interlayer interaction 
constants J", J'". This pJ|94el corresponds to the lattice 
structure of NH4CuCl3,t£HlZl with only nearest-neighbor 
exchange couplings taken into account. As before, we as- 
sume that all interactions are antiferromagnetic, and that 
J3 is much stronger than the other exchange constants, 
so that the system can be viewed as a set of weakly cou- 
pled dimers. Along the same lines of reasoning as those 
used in the previous Section, one obtains the effective 
spinless fermion model shown pictorially in Fig. ||; in ad- 
dition to (Q), there are now two more pairs of interaction 
and hopping constants, 

U" = - J", t" = 0, U'" = -J'" = -t'" . (13) 
2 4 K J 

There is again a particle-hole symmetry; the relation (|^) 
now has to be modified to 



One can find that up to the concentration M = | par- 
ticles can be pumped into the system without interact- 
ing with each other. The corresponding structure with 
M = 4 is shown in Fig. ||a; it is easy to see that no 
new particles can be inserted without interaction. Thus, 
the interaction energy Ei nt as a function of M always 
has a kink at M = \, and a symmetrical kink should be 
present at M = |. Those are exactly the values of M 
at which -the magnetization plateaus in NH4CUCI3 are 
observed.Ej However, it can be shown that in the atomic 
limit there are always other kinks at different values of 
M present in E- lrA {M ). Simply comparing the energies of 
different structures and displaying them in the i?i n t(M) 
plot, one can see that there always exists certain M = \ 
structure whose energy is below -Ei n t(^) + Ei ni (j), i.e., 
the corresponding point in the E[ nt (M) plot lies below 
the line connecting the M = \ and M — \ structures; 
in the same way one can check that M — i and M = | 
structures can be also favored under certain conditions 
on the model parameters. 

Since no other plateaus except M = | and M = | are 
observed in NH4CUCI3, this presents a problem. Let us 
consider a particular regime with the exchange constants 
satisfying the following inequalities: 



J" 1 < J", J' < 2J'", 
J' < 2Ji + J 2 + 4( J" - J" 



(15) 



/i «-> Mfcoie = -H -2U- 4*7' - 4(7" - AU" 



(14) 



Then one can show that in the atomic limit the only 
additional plateau structure realized in between M = 
1 and M = | is the M = \ structure shown in Fig. 

. Then, generally, there exist eight critical fields, as 
in the 2D problem considered in the previous section, 
which can be calculated perturbatively in the same way. 
The perturbation theory results in three dimensions are 
somewhat cumbersome even in the second order, so that 
they are listed in the Appendix. 

The important point is that in the regime determined 
by (|l^) hopping reduces the width of the M = i plateau 
already in the first order, while corrections to the width 
of M = j, j plateaus start from the second order. Thus, 
in the above regime there is only one additional plateau 
at M = i which is smeared by hopping more effectively 
than the plateaus at M — |, |. This could be a qualita- 
tive explanation of the real situation in NH4CUCI3, if one 
assumes that the additional plateau is completely wiped 
out. 



IV. DISCUSSION AND SUMMARY 

In this paper I show that in quasi-one-dimensional spin 
systems the location of plateaus in the dependence of 
magnetization on applied field can be considerably af- 
fected by presence of two- or three-dimensional inter- 
chain couplings. As an example, a simple model of two- 
dimensionally coupled zigzag chains is considered in a 
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dimer limit. It has been shown that, depending on the 
lattice structure and on the interplay between the model 
parameters, the system can exhibit plateaus at ^ and 
| of the saturation, in addition to the "usual" plateau 
at one half of the saturation determined by physics of a 
single chain. 

We argue that this mechanism can be rele- 
vant for the recept intriguing high-field magnetiza- 
tion measurements^ in NH4CUCI3, since the three- 
dimensional lattice structure of this material naturally 
favors the appearance of plateaus at \ and | of the satu- 
ration, i.e., exactly at those positions where the plateaus 
were actually observed. 

Unfortunately, the arguments presented here cannot 
be used for any quantitative predictions concerning the 
exchange constants in NH4CUCI3, because of the fol- 
lowing reasons: (i) zero-field magnetic ground state of 
NH4CUCI3 im plies that the first critical field H cl < 0, 
and from (A8) one can clearly see that then J3 cannot 
be much stronger than all the other couplings, as we as- 
sumed; (ii) to justify considering the hopping term as a 
perturbation, one would need to satisfy (|7|), and the con- 
dition J 2 — > 2 Ji does not seem plausible on the basis of 
the available crystal structure data.O 
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APPENDIX: GAPS IN 3D MODEL 



In the 3D model of Sect. [II the magnetization curve 
looks similarly to that presented in Fig. [| except that 
now we have plateaus at j and | instead ^ and |, re- 
spectively. The effect of hopping can be considered per- 
turbatively, provided that (0) is satisfied together with 
the additional inequalities 



J"<Ji, J"'<Ji. 



In the second order of perturbation theory, the dispersion 
s c z{k) of the excitation with one extra particle inserted 
above the M = | state, and the dispersion e C 2(k) of the 
excitation resulting from taking away one particle off the 
same M = \ state, can be represented in the following 
form 



s c<i (k) = e° ti + A, cos(fci + k 2 ) + B t cos 2(ki - k 2 ) (A2) 
+ 2d cos(fci — fc 2 )[cosfc3 + cos(fex + k 2 — ks)] , 

There are no first-order corrections to the energies of 
A/4 ± 1 states. The corresponding constants for the 
N / 4+1 excitation are given by the formulas 



t 2 



-c3 



H + 2U' + AU'" + — + — + 



2t 2 



U" U"> U' + 2U" 
6{t 2 + t' 2 ) 16[t' 2 + {t'") 2 } 



U -U' + 2U" - AU'" W + 3U'" 
I2t' 2 2t' 2 2t 2 



U> + AU" 



12f 



2U - W + 2U" U + U' + 2U" 
12t 2 At' 2 



U + 2U" U' + 2U" U + 3U" - U'" 

At 2 At 2 A(t"') 2 



(A3) 



-4, 



U' + 3U"-U"' U' + 2U"-U"' U' + U'" 

At' 2 2t' 2 16{t"') 2 

U + 2U" - U'" + U-U' + 2U" ~ 2U' + 3U'" ' 

m/_ A(t"') 2 

U -U' + 2U" - AU'" + U> + U'" ' 



(A4) 



B, = 



t'- 



C a = 



At'f' 



U" n ~° U' + W"' 
and for the A^/4 — 1 excitation one obtains 



8(t"') 2 , 2(t 

e c2 — h L ' 



t'- 



2t 2 



U' U" 
2t' 2 16[t' 2 + (t'") 2 ] 



2t^_ 

Z7 7 " + v + u" 

6t' 2 



U + U" 
12t' 2 



U' + 3U'" 
U(t'") 2 



U> + AU'" 2U' + 3U" 



U + 2U" 

6t 2 
TP + 2U" 



(A5) 



Ao 



Att' A(t"') 2 



f 



8t't'- 



U" 



W 



D SJ _ 

2 - jjm ' 2 - w + 3U „, 



Closing of the gaps A C] j = min£ Cj i(fc) (i = 2,3) deter- 
mines the critical fields H c2 , H c s marking the beginning 
and the end of the M — j plateau, respectively. Tak- 
ing into account that B 2 ^ and 62,3 are always positive 
[by our assumption on antifcrromagnetic character of all 
exchange interactions, see (ji)), (fl3"l)l, one can obtain the 
following formulas for the gaps: 

a_o / A + B-AC, B <C . 

-■■ + < A-B-2C 2 /B, B>C atA<u > 



(Al) and for A > the result is a bit more complicated, 



-A-B, 
a _ -n ) A + B-AC. 



C < {ABfl 2 
A<C,B <C 
B-A- 2C 2 /A, {ABfl 2 <C<A 
A-B-2C 2 /B, {ABfl 2 <C<B 



[We have omitted the subscripts 2, 3 above for the sake 
of clarity]. One can see that the expressions for A C 2, 
A C 3 "explode" when either of the interchain interaction 
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constants U', U", U'" vanishes; this corresponds to the 
obvious fact that the M = | plateau ceases to exist if 
any of the interchain interactions is absent. 

The energy e N/2+1 (k) = e { ^ /2+1 (k) + e { ^ /2+1 {k) of the 

excitation with one extra particle above the M = 4 state 
contains both first- and second-order corrections. In the 
first order approximation one obtains 



4/2+1 = P + 2 ( U + U ' + 2U ") 

+ 2t'cos(fci - k 2 ) + 2t"'[cosk 3 + cos(fci + fc 3 )], 
and the second-order contribution is 

(2 ) _ t 2 (l + cos2fci) 2t' 2 (l + cos2fc 2 ) 

£ N/2+1 - jj,, _ jjm + u -U' + 2{U" - U'") 

At 2 + At 12 + 8ttf cos(fci + k 2 ) 



(A6) 



U-U' + A(U" - U'") 
At 2 + At' 2 + 8tt' cos(fci - k 2 ) 
U + A(U" - U'") 
22t 2 



(A7) 



22t' ; 
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h u 
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At 2 


AU'" 4 
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f AU" - 
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u> - 


f AU" - 
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W" 


u + 


- W + A{U" - 


U'") + 


2U 


+ A(U" 


- U'") 



Zero of the corresponding gap A C 5 = xsxm.ej^i 2+ x{k), to- 
gether with the formulas (Q), [HI determines the critical 
field H c $ marking the end of the M = i plateau. 

The energy £ c i(k) of the one-particle excitation above 
the non-magnetic (M — 0) state is completely deter- 
mined by the hopping part of the Hamiltonian and has 
the form 

£ci(k) = /-!■ + 2tcosfci + 2t'[cosfc2 + cos(fci — ^2)] 

+ 2t'"[cosfc 3 + cos(fci + k 2 - k 3 )] , (A8) 

the gap in the above dispersion closes at the critical field 
H c i above which the magnetization starts to increase 
from zero. 

Exploiting the particle-hole symmetry (|lj) allows one 
to connect the following pairs of critical fields: 



Hs «-> H c i, H c5 <r-> H ci , 
H C 3 Hc6, H c2 <-> 7? C 7 



(A9) 



Setting H C 4 equal to H c5 , we can obtain a condition 
for "closing" the M = | plateau; restricting ourselves to 
the first order expressions for the sake of simplicity, we 
get 



Jl + 2 j2 



J" <J' + 3 J" 



(AlO) 



However, it is easy to see that ( AlO ) is already outside the 
range o f va lidity of the perturbation theory determined 
by (0), ( Al), so that this condition should be understood 
only as a rough estimate indicating that J', J'" have to 
be sufficiently strong to wipe out the plateau at M = h. 
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FIG. 1. Two-dimensional system of coupled zigzag chains 
described by the Hamiltonian (fe|) , a structure which i ij rnal ized 
in the (ac)-planes of the KCUCI3 family compoundsJi3t3 
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FIG. 2. Effective model of spinless fermions on a trian- 
gular lattice resulting from the system shown in Fig. |l|, as 
described by ([]). The structure with M = | in the "atomic 
limit" is shown; open and filled circles denote respectively 
empty and filled sites (i.e., dimers in |s) and \t+) states). 
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FIG. 5. Dimers coupled in a trellis lattice. This topology 
allows only one plateau at M = ^. 
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FIG. 3. Different processes of increasing the concentration 
in the atomic limit of the model ([|) at U' < U: (a) the 
process leading from M — J to the intermediate plateau at 



M = |; (b) the process increasing the magnetization above 
M = \. Half- filled and black circles represent respectively 
'new' particles inserted into the system and 'old' ones which 
formed the M = | or M — i structure. 
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FIG. 4. A schematic look of the interaction energy per 
site Ei nt £L function of the concentration M, in the "atomic 
limit" of (B). 




FIG. 6. Schematic plot of the magnetization M as a func- 
tion of the applied field H for the effective model (^). Solid 
and dashed lines show the magnetization curves for the cases 
U' >U and U' <U, respectively. 
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FIG. 7. Three-dim««si«nally coupled spin structure real- 
ized in KCuCl 3 familyHtJ Filled circles denote 5" = j spins 
of the lower layer in the (ac)-plane, and half-filled circles rep- 
resent spins of the upper layer. 
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FIG. 8. Visual representation of the effective spinless 
fermion model resulting from the model shown in Fig. ^| the 
interaction constants U and hopping amplitudes t are indi- 
cated for each link and are given by the formulae (§), @. 
Filled and empty circles represent occupied and empty sites; 
(a) the M — j structure; (b) the M = | structure which is 
realized under conditions (P-5|), 







